Fingerprinting Dark Energy III: distinctive marks of viscosity 



o 



0^ 
O 

u 

o 



> 

rn 

O 



X 



Domcnico Saponqj 
Instituto de Fisica Teorica (UAM/CSIC), Universidad Autonoma de Madrid, Cantoblanco 28049 Madrid, Spain 

Elisabetta Majerottqj 

Instituto de Fisica Teorica (UAM/CSIC) & Departamento de Fisica Teorica (UAM), 

Universidad Autonoma de Madrid, Cantoblanco 28049 Madrid, Spain 

(Dated: February 25, 2013) 

The characterisation of dark energy is one of the primary goals in cosmology especially now that 
many new experiments are being planned with the aim of reaching a high sensitivity on cosmological 
parameters. It is known that if we move away from the simple cosmological constant model then we 
need to consider perturbations in the dark energy fluid. This means that dark energy has two extra 
degrees of freedom: the sound speed c^ and the anisotropic stress a. If dark energy is inhomogenous 
at the scales of interest then the gravitational potentials are modified and the evolution of the dark 
matter perturbations is also directly affected. In this paper we add an anisotropic component to the 
dark energy perturbations. Following the idea introduced in [l[, we solve analytically the equations 
of perturbations in the dark sector, finding simple and accurate approximated solutions. We also 
find that the evolution of the density perturbations is governed by an effective sound speed which 
depends on both the sound speed and the anisotropic stress parameter. We then use these solutions 
to look at the impact of the dark energy perturbations on the matter power spectrum and on the 
Integrated Sachs- Wolfe effect in the Cosmic Microwave Background. 

PACS numbers: 98.80.-k; 95.36.+X 



I. INTRODUCTION 



The cause of the acceleration of the Universe's expan- 
sion, real or apparent it may be, is yet shrouded in mys- 
tery. The easiest explanation involves Einstein's cosmo- 
logical constant and results in the "standard" cosmolog- 
ical model, ACDM. However, this suffers of fine tuning 
problems, as well as all up to now available alternative 
models. Since the observations of supernovae type la 
(SNIa) 0, Qj many other probes have confirmed this ac- 
celeration, and several experiments arc planned in order 
to understand the nature of the phenomenon, usually 
called "dark energy" . The quality and quantity of the 
upcoming data will allow to better distinguish among dif- 
ferent models, both from the background expansion and 
from the cosmological perturbations point of view. In 
this context, it is important to determine and understand 
as well as possible all signatures characterising different 
models. 

In two previous papers [l|, |j] , some of us have studied 
the influence of dark energy pertturbations on the evolu- 
tion of cosmological perturbations and on the obscrvables 
related to it. To quantify this impact, the two parame- 
ters Q and S introduced in [5| were evaluated, the first 
measuring the fraction of dark energy density perturba- 
tions and the second related to the sum of the gravita- 
tional potentials, observable with weak lensing. In both 
[l|, |j] , dark energy was modeled as a fluid characterised 
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by an equation of state w = p/ p (where p and p are pres- 
sure and energy density, respectively) and a sound speed 
c^ = 6p/6p (where 5p are the pressure perturbations), 
both assumed to be constant. Paper [l| found a simple 
and useful analytical expression for matter and dark en- 
ergy perturbations, while paper [J] looked more carefully 
at observational implications. 

Here we take a further step in generalising the dark 
energy fluid. We introduce an additional degree of free- 
dom: the anisotropic stress. From the seminal work [y], 
we know that a classical scalar field has no anisotropic 
stress. Although most dark energy models are based on 
scalar fields (for a review see e.g. 0,01), it is nevertheless 
interesting, also based of the fact that very little is known 
about dark energy, to study this very general fluid and 
see how such a term would affect density perturbations 
and observables which are sensitive to them. Moreover, 
modified gravity models may be reformulated as effective 
dark energy fluids with anisotropic stress. The possibil- 
ity of detecting it is therefore related to the problem of 
distinguishing dark energy from modified gravity. 

Probably due to the absence of anisotropic stress in 
classical scalar fields, which are among the most popular 
candidates for dark energy, there is not much literature 
on this subject. One interesting paper on this is Q- Here 
a "generalised dark matter" with an anisotropic stress 
component was studied. This anisotropy corresponds, in 
the case of the fluid, to a viscosity term c'^^^ damping den- 
sity perturbations. The authors of this paper designed 
an equation governing the evolution of the anisotropic 
stress which recovers the free streaming equations of mo- 
tion for radiation up to the quadrupole. In [13] this same 
general paramcterisation of the dark energy component 



was analysed with data from the cosmic m.icrowavc back- 
ground radiation (CMB), large scale structure and super- 
novae type la, finding that both c^ and c^^^ were hard to 
constrain and that even future data would not improve 
very much in measuring them. A similar conclusion was 
reached by [ll| . Here forecasts were made on how well fu- 
ture CMB experiments will constrain an early, cold and 
stressed dark energy. Also 12[ used this paramctcrisa- 
tion to constrain extra neutrino species which are not 
explained by particles physics. A different approach on 
anisotropic stress was taken by |13| , which found general 
consistency relations for ACDM and studied anisotropics 
arising even in ACDM from second order perturbations. 
A non exhaustive list of other works concerning observ- 
able consequences of an anisotropic stress term is [ij, [l^ . 



In this context, our approach is similar as in [l| : we aim 
to try and solve analytically the perturbation equations 
for the simplest possible model of an anisotropic stress 
dark energy fluid in the simplest sensible approximations. 
This is useful both to understand more clearly the sig- 
natures induced by viscosity and to build simple tools 
for comparisons with observations. We parameterise the 
dark energy fluid by using three constants: the equation 
of state parameter w, the sound speed c^ and the viscos- 
ity cjjg. We start off from the anisotropic stress model 



of 0. 



The structure of the paper is the following. After dis- 
cussing the perturbation equations, defining our variables 
and describing the model for anisotropic stress in Sec. HIl 
we derive simple analytical solutions for the dark energy 
perturbations in presence of anisotropic stress during the 
matter domination era in Sec. IIIII We also verify that 
our expressions are a good fit to the numerical solutions 
obtained using CAME [la, [13] • In the same section we 
study the effect of the anisotropic stress of dark energy 
on the evolution of perturbations, concluding that the 
presence of a viscosity term c^j^, produces an anisotropic 
horizon, which results in an effective sound speed. In 
Sec. IIVI we compute the clustering parameters Q{k,t), 
S(A:, t) and /^(fc, t) defined in [^, for the case of our model. 
These have the advantage of tracking the numerical solu- 
tion even after matter domination, where they are strictly 
valid. Finally, in Sec. |V]we use our results to evaluate the 
effect of the viscosity c^,^ (and of the other dark energy 
parameters w and c^) on the amplitude and shape of the 
matter power spectrum, on the growth factor and on the 
integrated Sachs- Wolfe (ISW) effect. We stress that our 
aim is not to give a full theoretical analysis on imperfect 
fluid dark energy, but to analyse the effect of a possible 
viscous dark energy on the growth of perturbations. 



II. FIRST ORDER PERTURBATIONS IN DARK 
ENERGY 

A. Definitions 

In this paper we consider only spatially flat universes, 
and we use the Newtonian or longitudinal gauge so that 
our metric reads 



ds^ 



(1 + 2V') dr^ + (1 - 20) dxdx' 



(1) 



where a is the scale factor, r is the conformal time, if) and 
<f> are scalar metric perturbations (or potentials) and we 
do not consider vector and tensor perturbations. With 
H we indicate the Hubble parameter, H = {\/a){da/dt) 
computed with respect to the physical time t [dt ~ adr) 
and an overdot indicates a derivative with respect to r. 
Let us remind that while the gauge choice affects the 
perturbations on scales larger than the Hubble horizon, 
k<aH, on much smaller scales the observables arc inde- 
pendent of it. 

For a generic fluid with constant equation of state pa- 
rameter w = p/p, the perturbation equations are [iSl. Il9l| 



6' = 3{1 + w)(l)' 



Ha^ a 



Sp 



— wS 



(2) 



y - -(l-3^)^-F-g.^ + (l + H,^^+(3) 



Ha^ 



Ha^ 



(l + w) 



Ha'^ 



where S = Sp/p is the density contrast, V is the veloc- 
ity perturbation, dp is the pressure perturbation, a is 
the anisotropic stress of the fluid and with a prime we 
indicate a derivative with respect to a. Since we are in- 
terested in late times, we will assume that the Universe 
is filled only by two components: a matter fluid with 
w = dp ^ and a dark energy fluid, parameterised by a 
constant w ~^ — 1 , by an anisotropic stress function a and 
by a sound speed Cg[3J| related to pressure perturbation 
through 



Sp ~ clp5 + 



3aH lei 



^l) 



pV 



(4) 



where c^ — p/p is the adiabatic sound speed of the fluid. 
Here we restrict ourselves to models with w > — 1 to 
avoid crossing the phantom divide, where the above is not 
a good parametcrisation [20] ■ In our case the adiabatic 
sound speed is c^ = w since our dark energy fluid has 
constant w. 

Given the above assumptions the perturbation equa- 
tions ([2]) and (|4]) for our dark energy fluid become 



5' 



V 



V 



1 



9a2ij2 (c2 



fc2 



-{cl-w)S + 3{l + w) (/)' 



w)k'^ 



a^H 



a^H 



[^. 



(5) 
a]. (6) 



The first two Einstein equations (in conformal Newtonian 
gauge) are [18| : 



2 I J.' 



1 \ 47rGv-^ , 



where the sum runs over all types of matter that cluster. 
We can write Eqs. ([7]) and © in one single formula to 
obtain the Poisson equation 



i ^ 



ZaH 



v^ 



(9) 



where A^ is the (gauge-invariant) comoving density con- 
trast of the i*'^ fluid. The fourth Einstein equation, which 
describes the difference between the two potentials "0 and 
(/) as a function of the anisotropic stress, is 



fc2(0_V-) 



l2TTGa^{l + w)pa 
B{a)a 



^™o)a'''+'"Hl + w;)c 



(10) 



(11) 



where p is the energy density of dark energy. To be able 
to solve analytically the above equations, we restrict our- 
selves to the matter dominated epoch, which is also the 
most interesting regime in relation to observations. Here 
8TTGJ2iPi — 8TTGpm = 3iIo^mo«~^- A conscqueuce of 
this is that ^^ piAi ~ PmAm so that, at first approxima- 
tion, only matter perturbations contribute to source the 
gravitational potential in Eq. ^. Since we are interested 
precisely in the impact of a non-vanishing anisotropic 
stress of dark energy on the main observablcs, we should 
in principle not neglect a in Eq. (jlOp . However, we no- 
tice that the term B{a) in Eq. (|lip is proportional to 

1 + w ^ (with w 1) and to a^^^-^™) -- a"* which 

decays away very quickly. Hence, as a first approxima- 
tion, we neglect the contribution of a in the fourth Ein- 
stein equation, obtaining ip c:i (f>. 



B. Anisotropic stress daric energy 

To model a dark energy anisotropy is difficult as this 
might rise from a real internal degree of freedom of the 
fluid or from a modification of the geometry of spacetime. 
Scalar field dark energy for example always has cr = [g . 
However, while dealing with general fluids this might not 
be the case. We decide to look at the imperfect dark 
energy fluid model developed and analysed by Q and 
[lOl [21 |. In this fluid approach, the anisotropic stress 
is a viscosity term damping density perturbations. This 
anisotropy is sourced by scalar velocity perturbations and 



is generated by the shear term in the metric fluctuation; 
moreover it has to be gauge invariant ^ . 

In this class of models the dark energy anisotropic 
stress cr is a function of scale k and time a. The equation 
for a is given, in Newtonian gauge (see (Sl), by 



V 



3{l + wya^H 



(12) 



where the viscosity parameter cf,^^ quantifies the coupling 
of the anisotropic stress to the velocity perturbation. On 
the left hand side of Eq. (fT2|) also appears a Hubble 
drag term. This equation, governing the evolution of the 
anisotropic stress, is built to recover the free streaming 
equations of motion for radiation up to the quadrupole. 
We remark that this ansatz does not cover mainly models 
in which stress fluctuations are not sourced by density 
and velocity perturbations but act as external sources 
for the perturbations J22| : here also vector and tensor 
stresses have to be modelled in order to understand the 
evolution of perturbations. 



III. SOLUTIONS FOR THE DARK ENERGY 
PERTURBATIONS 

If the universe is matter dominated, then k'^ip ~ 
—4:TrGa^PmAm — — 3/2 TJgrimo'^o = ~</'o is constant and 
the dark matter density contrastgrows linearly with the 
scale factor. Along the lines of [l[, we look for solutions 
for S and V in presence of an anisotropic stress evolving 
as in Ea.(|12p in two regimes: that of perturbations below 
and of those above the sound horizon. 



A. Modes below sound horizon 

We start with perturbations below the sound horizon, 
i.e. such that k <C aH /cg. Here, we can neglect the terms 
proportional to V and V appearing in Eq. ^ as they 
arc small compared to the third term on the right hand 
side, which is ^ fc^c^- We are left with an equation for 5: 



6^ 



1 



fc^c? 



(00 + fc^cr) 



(13) 



where we used k (j) ~ — 0o- Substituting Eq. (|13p into 
Eq. ([5]) we obtain another differential equation for cr, cou- 
pled to V: 



V 



{cl -) 



1 



w 



fcV) (14) 



where we have neglected the second term in square brack- 
ets of Eq. ([S]), as fc~^ is small for modes below sound 
horizon. Inserting the expression for V of Eq. p4l) into 
Eq. (J12p . we obtain a first order differential equation in 




FIG. 1: Behaviour of the dark energy density contrast S for 
scales above sound horizon as a function of the scale factor a. 
The black dot-dashed line represents the numerical solution 



computed with CAMB for a model with c„i 



10" 



and w = —0.8 for the mode k = 200-ffo (where radiation was 
omitted in order to obtain a longer interval in a and show 
all different dynamical regimes). The red solid line is the 
approximated solution for c^^g — for scales above sound 
horizon, Eq. (|27|l . while the blue dashed line is Eq. (|26p . valid 
for c^ig 7^ and scales above sound horizon. The vertical 
dotted lines correspond to the value of a at which the mode 
enters the causal horizon (left brown line) and the anisotropic 
horizon (right green line) . We see that the numerical solution 
consistently follows first the red solid line and then the blue 
dashed line once the anisotropic horizon is crossed. At a ~ 1 
we start to see a deviation from our approximated solutions 
due to the end of matter domination, hence of their domain 
of validity. 



a only 



1 31 + u; 






3c2 



w 



(15) 



a h?c1 
This equation can be solved exactly, the solution being 



"^vis V^s 



3c2(l + H+8(c2-^)cLfc^' 



(16) 



plus an extra term rapidly decaying when a grows (see 
Appendix |X|. Substituting a into Eqs. P^ and P^ we 
finally obtain expressions for S and V: 



3{i + wy 



V = - 



3c2(l + u;)+8(c2-u;)cLfc' 
9(1 + wf (c2 - w) 



-3aH {cl -w)S. 



(17) 



(18) 



We notice that in our approximations a and 5 arc ap- 
proximately constant while V ^ a~^'^. 



Dark energy perturbations in this scenario are sup- 
pressed by a factor (1 -I- w)^. The (1 + uif' factor is 
originated both by the {\ + w) factor contained in the 
gravitational potentials, which source dark energy per- 
turbations, and by the anisotropic stress. Therefore, 
anisotropic stress slows the growth of dark energy per- 
turbations by an extra factor of (1 -h w) with respect to 
the standard case of clustering quintessence. In the limit 
of no viscosity, c^j^ = 0, we consistently obtain a ~ 
and 



5 = (1 + u;)- 



,21,2 



V 



2.{l + w){cl-w)HoyJV[^, 



-1/2 



s^fc-^ 



(19) 

(20) 



which corresponds to the standard case studied in [l|. 
Moreover we see that now the suppression factor of dark 
energy perturbations is again (1 -I- w) and not anymore 
{l + wf. 



B. Modes above sound horizon 

To find solutions for a^ 6 and V for modes above sound 
horizon, we must again start from the velocity pertur- 
bation equation ([6]) but this time we neglect the term 
proportional to c^ so we have 



v 



-V = — J-— (00 



fcV 



(21) 



This equation, together with Eq. ([T2|) . calculated in 
the approximation of matter domination i.e. for a?H = 
Hq yj^mo V^-t forms a closed system that can be solved 
after some mathematical manipulation, giving the follow- 
ing solution for a 



1 



3 

2aa 



9 



2(aa)2 4(aa)3 



where 



8cl 



a 



O-tlQ i Ljjl 1 



(22) 



(23) 



plus some extra decaying terms which we can neglect. As 
in the case of modes below sound horizon, the dominant 
term is a constant. Looking again at Eq. (pTj) . we now 
notice that, contrary to the standard case where c^j^ = 0, 
here too we could have neglected the terms ex V,V', since 
they decay as a~^/^ compared to a. Indeed neglecting 
them, Eq. ([21]) gives us consistently a = —(pg/k^. This 
also coincides with the solution for a we found for modes 
below sound horizon (where we also had neglected terms 
proportional to V, V), if we set c^ = into our solution, 
and the same happens to V , which is 



V = - 



9(1 + u;)^ 
8c^- 



-ffov^^ 



'P 



-1/2 



(24) 



We compute S from Eq. ([5]), where we neglect terms 
and ~ k~^: 



= —^5 — TT 

a 8c7,- a fc^ 



and the (dominant term in the) solution is 

8ci„w F ' 



(25) 



(26) 



It is clear again that the last equation is a subcase of 
solution ([TT]). Looking at Eqs. (|22j|23p it would seem that 
the terms containing c^j^. may come to dominate if we 
are dealing with small values of c^j^, especially at early 
times. However, the impact of the these terms on the 
growth of the dark energy density contrast is negligible, 
as explained below. 

Let us remind that these results have been obtained 
under the assumption of a time- independent w, c^ and 
c^ j3 . but we have not excluded a fc-depcndence of c^ and 
'^vis- ^^ ^^ ^^^^ important to notice that the solutions 
found above do not reduce to the case c^ = c^^ == 0. For 
this case we need to use the solution found in [l| for the 
density perturbation: 



(5 = (5o(l + w) 



0/7q i l^Yif.^ 



l-3w 



(27) 



In Fig. [T] we plot the numerical solution for the dark 
energy density contrast for k = 2007Jo as well as the an- 
alytic solutions for modes above sound horizon (c^ = 0), 
Eqs. dMl) and ^ for cl,, = lO"" and c^-^ = respec- 



tively; here and in the following figures we set fir, 



0.23 



and Hq = 70 km/s/Mpc. The values are chosen so 
to show the full complexity of the 5 evolution. It can 
clearly be seen how perturbations start to grow as soon 
as they enter the causal horizon. Although in this case 
there is no sound horizon as the sound speed is zero, the 
presence of a non-vanishing anisotropic stress creates a 
new anisotropic horizon in the evolution of the dark en- 
ergy perturbations: when perturbations cross this hori- 
zon then the viscosity of the fluid counteracts the grav- 
itational collapse and prevents the perturbations from 
growing. The value of the scale factor at which pertur- 
bations cross the anisotropic horizon can be evaluated 
from kcyis ~ aff|35l|. 

It is the presence of this anisotropic horizon that jus- 
tifies the assumptions made so far. To illustrate this, 
let us consider again the solution for a obtained assum- 
ing Cg = 0, i.e. Eq ([2^ . If the sound speed is zero 
then the only component capable to prevent dark en- 
ergy perturbations from growing is the anisotropic stress, 
which starts becoming important and counterbalancing 
the gravitational collapse only after the anisotropic hori- 
zon has been crossed by 6. A way to understand this 
is by looking at Eq. P^ : when c^^^ is very small, the 
evolution of a is essentially decoupled from V, hence the 
term proportional to a in Eq. (j2ip is negligible and the 



anisotropic stress does not affect the growth of perturba- 
tions. Only later in time, approximately when crossing 
the anisotropic horizon, does viscosity enter the game. 

As shown, the size of this horizon depends on the 
value of the viscosity term c^.^g. During matter domi- 
nation, the anisotropic horizon crossing happens at a ^ 
H^^mo / {k'^ c'iis)'- the smaller is c^j^, the later in time 
perturbations enter the horizon and consequently they 
have more time to grow. The reader might think that the 
terms appearing in Eq. (|16p come to dominate eventually 
for a small value of c^^^, when a is small. However, this 
is not the case thanks to the presence of the anisotropic 
horizon. This is illustrated in Fig. [T] where c^ = and 
c^jg assumes a very small value, c^j^ = 10"'* . Before en- 
tering the anisotropic horizon dark energy perturbations 
grow as if the anisotropic stress were absent, following 
Eq. (|27p . Only when the anisotropic horizon has been 
crossed perturbations stop growing because a becomes 
important. However, this only happens at late times, as 
c^j^ is small. Hence, being a large, only the first con- 



stant term in Eq. 



is important. To prove this, let us 



consider the decaying mode in Eq. 
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10" 



10- 



(28) 



where k ~ 200i/o- The above term is larger than the 
constant term if a < 10"^, but it can never dominate 
because the anisotropic horizon is at a ~^ 10^*. 

What happens at large scales when c^^^ is larger, say 
c^j3 ~ 10~^? From Eq. (|22|) it would seem again that 
a cannot be approximated by a constant. Here a differ- 
ent mechanism comes into play: large scales enter the 
causal horizon later than small scales (e.g. during mat- 
ter domination at a ^ H'^Vlmg/k'^), hence perturbations 
with small k cross the causal horizon at sufficiently large 
a so that the time-dependent terms in Eq. p2|) can be 
neglected. If e.g. we take scales as large as fc = Hq, the 
causal horizon is at a = ^ma ^^^^d here we get (again for 
w = -0.8) 



1 
aa 



3(1 



0.2 < 1. 



(29) 



This quantity becomes even smaller for values of w closer 
to — 1 . The considerations made in this section ensure us 
that the time dependent terms (pS)) decay fast enough in 
all relevant cases. 



C. An effective sound horizon 

It is interesting to notice that the anisotropic stress is 
not always the dominant component affecting the growth 
of dark energy perturbations. As an example, we plot in 
Fig. ([2|) the evolution of the dark energy density con- 
trast 5 for different values of the viscosity term c^^^ 
(with the corresponding anisotropic horizon) while keep- 
ing the sound speed cj. fixed. As the viscosity term dc- 



creases, perturbations start to feel the presence of pres- 
sure perturbations, which, similarly to c^^^, damp S once 
it crosses the sound horizon. At this point, having density 
perturbations already been damped, the anisotropic hori- 
zon plays no role anymore and the viscosity term becomes 
unimportant. This effect can be clearly seen in Fig. ([2|, 
where the red solid line and the green long-dash-dotted 
line, corresponding to c^^^ = 10~^ and c^^^ = lO""*, re- 
spectively, are basically unaltered when crossing their 
anisotropic horizon (that for the value used it happens at 
a ^ 10~^), having already experienced the damping due 
to pressure perturbations. A similar behaviour can be 
seen when the anisotropic horizon is crossed before the 
sound horizon: here S is damped by the viscosity of the 
fluid and the crossing of the sound horizon later in time 
produces no effect. Furthermore, looking at Eqs. pTIIlSp 
we notice that it is possible to rewrite them in terms of 
an effective sound speed: 



^eff 



3 (1 + w) ' 



(30) 



We list the equations for 6 and V in terms of c^jf and 
Cg in Appendix O In view of these considerations, the 
degeneracy between cj and c^^^ found by us, [l^ and 
[ll| is confirmed and has a clear explanation. From c^fr 
it is hence possible to define an effective sound hori- 
zon, which takes into account both c^j^, and c^ and de- 
termines the tiine at which perturbations are damped 
by the anisotropic stress or by pressure perturbations: 
a?H^ ~ A;2[c2 + ^/2,cl,,{cl - w)l{l + w)]. We can ask 
when the anisotropic stress comes to dominate over pres- 
sure perturbations. The sound speed starts becoming 
more important than anisotropic stress if c^^^ < c^/10 
(for our value of w), which is also confirmed by numeri- 
cal results. 



IV. CLUSTERING PARAMETERS 

When dark energy starts dominating, starts decay- 
ing; moreover, dark energy is also clustering. This vari- 
ation can be conveniently represented using the param- 
etcrisation defined in |5| where two new quantities were 
introduced: Q, which is directly related to the Poisson 
equation and describes the amount of clustering dark en- 
ergy, and 77, which is related to the amount of anisotropic 
stress. Using this new parameterisation, the Poisson 
equation reads 



fc20 = -Ai^Ga^QpmX, 



where 



g-1 






(32) 
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FIG. 2: Numerical evolution of the dark energy density per- 
turbation 5 with k — 2QQHij for different values of c^j^. The 
black dotted, the brown long-dashed, the blue dot-dashed, 
the green long-dash-dotted and the red solid line (bottom to 
top) correspond to cl^^ = 10"\ 10"^, 10"^, 10"* and 10"^, 
respectively. The vertical lines give the scale factor at which 
each mode of matching colour/line style enters the anisotropic 



horizon (left to right), 
was assumed. 



For all solutions, w = —0.8, c„ = 10 



measures the amount of clustering dark energy whereas 
Q refers to the sum of dark matter and dark energy clus- 
tering. The parameter 77 is instead defined from 



4) = {l + Tj{a,k))(j>. 



(33) 



We can now compute the gauge-invariant density per- 
turbation A of matter and dark energy. During matter 
domination we have at first approximation A„i = Jq a, 
where 5q is a constant of integration. For dark energy we 
obtain 



(1 + w) 00 



(34) 



■-e// 



where we have neglected the term proportional to V since 
it is decaying. Also recall that this solution is valid for 
modes both above and below sound horizon. 

In our case, we will find an expression for Q — 1 valid 
for aU c^.^ > (and c^ > 0): 



1 



3gg (1 - r?„, J (1 + u;)a- 
2/c2 



(35) 



•-eff 



We notice that the above equation is not valid for the 
C3X) special case c^ = c^j^ = (as it was already the case for 
5) . Being the latter the standard case with no anisotropic 
stress, it had already been calculated in [l[ to be 



Q{cl 



0)-l 
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1 -n 



1 — iw 51^ 



"'" n^^^" 



(36) 



The reason why we mostly express quantities in terms of 
Q — \ rather than Q itself is because the former directly 



To obtain a single expression for Q, including anisotropic 
stress as well as the case c^ = c?„„ = and valid for modes 





FIG. 3: Behaviour of Qtot — ^ for two difTerent values of the vis- 
cosity term c^^^. The sohd hues are numerical solutions com- 



FIG. 4: Behaviour of 
the viscosity term c^^j. 



"^tot — 1 for two different values of 
The solid lines are CAMB outputs 



puted with CAMB (blue/upper line: c^is = 10 , black/lower (blue/upper line: c^^^ — 10 , black/lower line: c^^^ — 10 



line: 



= 10" 



while the red dashed line and the yel- while the red dashed line and the yellow dot-dashed line repre- 



low dot-dashed line represent Eq. (|37|l with c^^^ = 10 * and 
Cvis = lO"'^, respectively. All lines correspond to fc = 200i/o, 



w — —0.8 and c. 



10" 



sent Eq. (|37p with c^^^ — 
All lines correspond to k 



10 and c^„ 
= 200Ho, w ■ 
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10" 



above and below sound horizon, we define 



[Q-i][Qic 



0)-l 



i-a. 



iin 



[Q-i] + [Q{c 

'-{l + w) 



2 _ 
vis 
. — 5w 



0)-l] 



l-3w 



2k^c 



(37) 
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The above expression is extremely simple and it looks 
hke the solution found in [l| , except for the term includ- 
ing the effective sound speed c^ff which depends both on 
the sound speed and on the viscosity term. It is easy to 
see that if the anisotropic stress is absent then Eq. ([57)1 
reduces to the expression found in [l[. From Eq. ([37)) 
we have a further confirmation of the degeneracy of c^ 
and c^j^, when looking at quantities which depend only 
on the clustering of dark energy. To break this degen- 
eracy, we will have to combine measurements of Q with 
measurements of 77, defined in Eq. (j33p . as will be clear 
from the following. 

Let us indeed evaluate the quantity 77: 
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where we used the relation k 
rameter depends not only on 
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'-eff , 



(38) 



e// 



-4'aQtot- This pa- 
but on the ratio of 



the sound speed and c^fr- For this reason, the only way 
of breaking the degeneracy between c^j^ and cf (at least 
at the linear perturbation level) is through the combina- 
tion of observables constraining Q and ry. The detection 
of a non-zero rj alone would prove the existence of an 



anisotropy but it would not allow to quantify it as only 



the ratio 



ll(?eff would be measurable. 



If instead of 77 we prefer to use the variable E intro- 
duced in [5| and more strictly related to the weak lensing 
potential ^I* = + 0, which is defined as 



E = 



1+2^ 



(39) 



we can evaluate it analytically using Eqs. ([57)1 and ([55]): 
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(40) 



where /3 = 2k'^clf^/[iHlVLr,^„{l ~ 3u;)]. The above equa- 
tion shows that the weak lensing parameter E differs from 
Q by about 10 — 15% (depending on a) in case of a non- 
zero anisotropic stress. 

In Figs.[3]and|3]we compare Qtot and Titot, respectively, 
with their numerical solution, computed using CAMB, 
for different values of the dark energy viscosity term keep- 
ing the sound speed fixed. We use two different values of 
c^jg that are respectively ten times smaller and ten times 
larger than c^ ; these values were chosen such that in the 
first case the sound speed is the dominant component 
while in the second the viscosity term dominates. We 
find very good agreement between the numerical solution 
and our analytical estimate, as it was the case in [l|. 



V. EFFECT ON SOME OBSERVABLES 

A. Matter power spectrum 

As in [l[, after evaluating the impact of anisotropic 
stress on dark energy perturbations, we move to look at 
the (much smaller) effect it should have on matter pertur- 
bations. These are affected by the change in gravitational 
potential from (pm = — 3/2(5oiifo^mo/^^ to (j) ~ 4>m + 4''DE, 
where 0de = (f'miQ ^1)- For modes larger than the 
dark energy sound horizon, Q is given by Eq. (j35]) where 
c^ = and we can solve for the dark matter velocity 
perturbation, Eq. ^ with 5p„i = 0, to find 



Vm = —Hoy'^lmo^OVa 
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-3tu-l 



1/3 -2i 



where 






(1 



(41) 



(42) 



and where the decaying term has been omitted. This 
equation is of course strictly valid only during matter 
domination, but it gives again the right order of magni- 
tude of the effect even after this epoch, as we will see 
in the following. The correction factor with respect to 
the case of no dark energy perturbations is the term in 
square brackets. We notice that its evolution in time, 
contained in the second term in square brackets and pro- 
portional to a~^'^~^, is different from that in absence of 
anisotropic stress, which was r^ a~^"'. Substituting Eq. 
(PT|) into Eq. ^ (clearly with w,„ = 6p,n = 0) we obtain 
a differential equation for (5,„, whose solution is 
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(43) 



As in the case analysed by [l|, the enhancement factor 
due to the viscosity term varies depending whether per- 
turbations arc inside or outside the causal horizon. In the 
first case and for perturbations outside the sound hori- 
zon, the first term in curly brackets is more important 
so that the correction term is Qoa~^'^ /[w{6w — 1)], while 
for perturbations outside the causal horizon, the term is 
9H^n,n„Qoa-'''"-^/i2k^). 

Let us now compare the matter power spectrum 
P{k,a) ^ 6m{k,a) obtained in the case of c^jj, 7^ to 
the standard one, obtained in absence of dark energy 
perturbations P^'^°(fc) ~^ Sf^^{k,a). In Fig. [5]we plot 
the ratio of these two quantities. As can be seen, our 
approximation of dm is of the same order of magnitude 
of the numerical solution, which is a good result for a 
second order quantity such as Sm- We also see that once 
the matter perturbations enter the anisotropic horizon, 
the analytical and numerical solution coincide and both 
tend to the solution of unclustered dark energy. This 





1.05 




CO 


1.04 


s I : 






N 1 : 


M 


1.03 


\ ■ : 


Q 




\ : : 


m e 


1.02 


\ 1 : 


M 




^^^k 1 J 


^ — ^ 


1.01 


^^^L 


h 




^T^k 


^ 




■^^^^fc 




1.00 


^^^^fc_^^ 








0.99 


I i 



0.01 



0.1 



k[h Mpc ^] 



FIG. 5: Behaviour with k of the ratio of 5^, to the square of 
the standard matter perturbation in absence of dark energy 
perturbations [5^^°]^. The black solid curve represents the 
numerical solution found using CAMB with c^^^ = 5 x 10~^ 
and C3 = 10~®. The blue dashed curve is our analytical solu- 
tion, computed from Eq. (|43p . The green dot-dashed vertical 
line represents the anisotropic horizon while the red dotted 
line is the sound horizon. Here we fix ui = —0.8 and a — 1. 



is consistent with what we have learned from the previ- 
ous section, i.e. that dark energy perturbations inside 
the anisotropic (or effective sound) horizon are damped 
and can be neglected. An implication of this is that 
anisotropic stress will be difficult to detect being present 
only at large scales. Also, for w = —0.8 we expect a 
maximum 4% enhancement of P(fc) on scales larger than 
both the anisotropic and the sound horizon. Our numer- 
ical calculation shows that the enhancement is closer to 
2%, hence the effect is very small. 

Having evaluated the impact of the anisotropy on the 
matter power spectrum amplitude, we move to compute 
explicitly a more easily observable parameter. 



B. The growth factor 

In the ACDM scenario, dark matter perturbations 
grow logarithmically with the scale factor a during ra- 
diation domination and when matter domination starts 
they grow linearly; at late times, when dark energy starts 
to dominate, the growth of dark matter perturbations is 
suppressed. It is known that in the ACDM model the 
growth factor can be defined as: 
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(44) 



where 7 ~ 0.545 is called the growth index. There are 
two ways to modify the growth factor. First, one can use 
a different background cosmology, associated to a differ- 
ent Hubble expansion. Second, perturbations can differ: 
if dark energy starts to cluster, then the gravitational po- 
tentials change, being sourced by the total density per- 



turbations and by the anisotropic stress (see the Poisson 
equation © and the fourth Einstein equation ([TU|) and 
this will affect the growth rate of dark matter. It is pos- 
sible to include all these effects in the growth index 7 
and we therefore expect 7 to be a function of w, c^ and 
clig{oT equivalently of w, Q and 77). 

In [2J| , the dependence of the growth index from the 
dark energy clustering parameters Q and 77 was described 
as 
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(45) 
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where 



A{Q,v) = 



il + v)Q-l 
1 - fl,n (a) 



(46) 



While dark energy perturbations themselves are difhcult 
to measure unless dark energy has very low sound speed 
and viscosity, the growth index seems to be a more easily 
detectable parameter and several ongoing and future ex- 
periments are built to measure its value. It was shown in 
[J] that the presence of a sound speed (when c^j^ = ry = 0) 
in the dark energy perturbations always decreases the 
value of the growth index 7 because Q — 1 is always posi- 
tive due to the relative increase of dark energy perturba- 
tions. Here instead the addition of a non-zero anisotropic 
stress may lead to an enhancement of the growth index. 
To prove this let us consider a particular case. If we 
take modes above the sound horizon, using Eqs. (j35]) and 
Eq. (p8|) with the sound speed set to zero, Eq. (|46l) reads 
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The fmiction B{a) is always positive as long as if is larger 
than —1, hence the term in brackets can be negative only 



if c^jj > (1 -I- w)/{—8'w). If we assume 



-0.8 then 



A{Q, 77) is negative as long as c^j^ > 1/32 and 7 is en- 
hanced. However, the effect of the anisotropic stress on 
scales of interest is very small due the smallness of the fac- 
tor B{a)/k^ in Eq. (PS|) that contains the term {l+w)/k'^. 
Again, if the anisotropic stress is zero, then dark 
energy perturbations always decrease the value of the 
growth index; we find that the strongest deviation from 
IKCDM = 0.545 happens when c^ — and c^j^ ~ and 
it is of the order of 3%, see [l| . If we want to increase the 
growth index and at the same time we also want pertur- 
bations in the dark energy sector then we need a non zero 
anisotropic contribution; this is the case for the DGP 
model (treated as an effective anisotropic dark energy, 
see [2J]) where the growth index is 7 ^ 0.68. However, 
the type of anisotropic stress used in this paper is un- 
able to increase the growth rate 7 from jacdm to such a 
large value. For instance, if we assume the viscosity term 
to be c^jj = 1 then A{Q,t]) ~ —1.5 x 10~^ for scales 
k ~ 200iJo- 



FIG. 6: ISW power spectrum for cj^s = (red solid line) and 
Cvis ~ 10~^(blue dashed line) fixing w — —0.8 and c^ = 10~*. 



C. The Integrated Sachs- Wolfe effect 

Another potentially observable effect of anisotropics 
is the integrated Sachs- Wolfe (ISW) effect. Indeed, as 
explained e.g. in j25H28l |. being this a late-time effect, it 
is the part of the CMB spectrum where dark energy has 
the largest impact. The total temperature shift of the 
CMB photons due to the ISW is ^■. 
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To 

aH-dx 
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dr 



(48) 



where x is the comoving distance {dx = —cdr = —cdt/a). 
In this work we decide to consider the effect of the 
anisotropic stress only for dark energy, in order to iso- 
late its contribution and to compute its size analytically. 
See llj for a numerical analysis which includes massive 
neutrinos and their anisotropic stress. The term inside 
the integral is the derivative of the weak lensing poten- 
tial ^ = tl) + <f> with respect to the scale factor, which in 
Fourier space is 



$' 



|^g^{S(a,fc)A;,(a,fc) + 
1, 



+ S' (a, k) A,„ (a, k) - -S (a, k) A„ (a, k) ^ . (49) 

From this expression we can see that anisotropic dark 
energy perturbations modify the ISW effect through the 
changes induced in A,„ and through the additional pres- 
ence of E and S'. 

At linear order, it is possible to isolate today's A^ 
from its time evolution: 



A„ (a, fc) = aG (a, k) A,„,o (k) , 



(50) 



where A,„ (k) = A„i(a = 1, fc). Already in jl| and [301 ^ 
deviation from the usual scale- independent growth factor 
G{k) due to dark energy perturbations was evidenced. 
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There, the fc-dependence was determined by the presence 
of the sound horizon. In our case, the mechanism is the 
same, but depends on the effective sound horizon, which 
is determined by both the sound speed and the viscosity 
of dark energy. Here we express in first approximation 
G{a,k) = an(a)''("'''), with 7(a,fc) from Eqs. (|45|) and 
(|4S| . We write Eq. (09]) as: 



$' 



3 Hq r2,„o d 



2 fc2 da 
Eq. (|i5|) reads now 



{G(a,/c)S(a,fc)}A„,o(A:) . (51) 
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(52) 



3 Hn^r. 
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^a^iJ^JG (a, fc)E (a, fc)}. (53) 



Hence the ISW-auto correlation spectrum Cc^c, (^) is (in 
the Limber-projection [3l| and in the flat-sky approxi- 
mation) : 



Ca (e) 



'^^ , wHx) 
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(54) 



where Paa {k) is the hnear matter power spectrum to- 
day: 



fc^-PAA (fc) _ c2 f k 
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T2 (k) , 



(55) 



6h is the amphtude of the present-day density fluctua- 
tions at the Hubble scale and T (k) is the dark matter 
transfer function. 

In Fig. m we plot £{£ + l)Ce/{27r) using Eq. ^ for 
two different values of the dark energy viscosity keep- 
ing the sound speed fixed to 10""'. For T{k) we use the 
fitting formula of |32| . which agrees with the numerical 
results from CAMB with a precision sufficient for our 
purposes. In computing the matter power spectrum we 
neglect the effect of dark energy perturbations evaluated 
in Sec. IV Al having ascertained that it is very small. This 
means that the two curves corresponding to different c^^^ 
differ because of the term {d{Y,G) / da)'^ . 

Let us then evaluate 



(GS)' ^G'S + GE'. 



(56) 



While both S and S' differ from the case with no pertur- 
bations, corresponding to S = 1, S' = 0, the first devi- 
ates too little from 1 to explain the differences between 
the red and the blue curves of Fig. |6l We have checked 
this by fixing S = 1 and letting E' 7^ and noticing that 
the results are unchanged. This means that the effect of 
c^^g on the ISW comes mainly from the term E'G. Using 
Eq. (gni) and ^ we find 

S' = g' + ^B'{a) (X - (S - 1) (57) 

200 a 
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FIG. 7: Magnification factor A of Eq 



ent values of the viscosity term: c^^ 
and c^ig = 10~^ (blue daslied line). 



for two differ- 
— 10^* (red solid line) 
Here we fix the scale of 
perturbations to fc = 200-H'o and our dark energy model has 
-0.8 and cl = 10"*. The vertical lines at a = 0.081 
0.24) show the time at which the perturbation with 
^ {ciis = 10^*) enters the anisotropic horizon. 



w 
(a 
c^ =10 



where the last proportionality holds because (Q — 1)' oc 
{Q - 1) /a and B'{a) oc B{a)/a hence aS' ex (S- 1). Let 
us now go back to Eq. ((56)) and look at the final total 
effect. We know that as dark energy slows down the 
growth of dark matter perturbations, at late times G' is 
always negative. E' instead can be positive or negative 
depending on the value of c^j^ , although in both cases it 
is a small number. Therefore we have two possibilities: 
the first is when the two contributions partially cancel 
and so make the effect smaller (this is in general the case 
for small values of the viscosity term); the second is when 
the two contributions sum up enhancing the total effect. 

To visualise this we proceed as in [l| by defining a 
magnification parameter for the ISW power spectrum 



A' 



d{G{a,k)T.{a,k)) /da 



(S 



dC (a) /da 

Y.'G/G'f, 



(58) 
(59) 



which we plot in Fig. [7] for two values of c^^^ for k = 
200i^o- Here we keep the sound speed and the viscos- 
ity term small enough so that the dark energy pertur- 
bations have time to grow before entering the effective 
sound horizons. For small c^j^ the dark energy perturba- 
tions are able to grow for sufficiently long time and they 
partially cancel the contribution from G' and decrease A 
by about 40%. As we increase the viscosity term c^^^, the 
effects become smaller because the dark energy pertur- 
bations cross the effective sound horizon early, increasing 
A. Based on these considerations, we expect again the 
viscosity term to act similarly as the sound speed: a low 
value of cjjg decreases the power of the ISW effect, dif- 
ferentiating it from the standard case. 

As mentioned before, there might be cases where E — 1 
is negative, enhancing the total ISW effect. In Fig. ([5]) 
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FIG. 8: Magnification factor A of Eq. (l58|) for two differ- 
ent values of tfie viscosity term: c^^^ — (red solid line) and 
Cvis = 1 (blue dashed line). Here we fix the scale of perturba- 
tions to fc = 2OO//0 and our dark energy model has w = —0.8 
and c? = 1 



wc show the magnification A for cf,^^ — 1 and Cy^^ = 
0, setting c^ = 1 in both cases. If the sound speed is 
large enough then S' is negative and it sums up to G" in 
Eq. ([56|) . Nevertheless the effect is too small (of the order 
of 10^^) and it is practically unobservable in the C^s. 



VI. CONCLUSIONS 

In this paper we have studied an imperfect fluid dark 
energy with no n- vanishing viscous anisotropic stress. 
The model we have considered is described by its back- 
ground constant equation of state u; < and at lin- 
ear order in perturbations by its sound speed c^ > 
and its viscosity parameter c^^^ > 0. For the dark en- 
ergy density and velocity perturbations we have found 
simple analytical expressions represented in Eqs. (J17p 
to ((26| . Although we have assumed matter domination 
throughout this work, the functions Q{k,a;w,cl,cf,^g), 
'i'l{k,a;w,cl,cl^g) and I](fc, a; w, c^, c,^,,,) which can be 
computed with our simplified analytical expressions are 
approximating well numerical results even at later times, 
when dark energy starts dominating. Our analytical 
expressions can be seen as the "fingerprints" of a vis- 
cous dark energy model (or of a modified gravity model 
that can be effectively described by such an anisotropic 
stress) . 

In particular we find that the dark energy perturba- 
tions are always smaller than the perturbations in dark 
matter, by a factor (1 + w) if the viscosity term is small, 
'^vis '^ ^ ^^'^ ^y ^^ even larger factor of (1 + if )^ when 
the anisotropic stress is not negligible. 

If dark energy has zero anisotropic stress then it is pos- 
sible to divide the evolution of its perturbations in three 
phases: outside the causal horizon, inside the causal hori- 
zon but outside the sound horizon and inside the sound 
horizon. In the latter regime their growth is suppressed. 



In the presence of a non zero viscosity term, a further 
horizon comes into the play, which we dubbed anisotropic 
horizon. Once perturbations enter this horizon, they are 
damped as in the case of the sound horizon. This adds 
in principle two extra regimes: perturbations larger than 
the anisotropic horizon and perturbations smaller than 
the anisotropic horizon. The last two regimes may hap- 
pen before or after the sound horizon depending if c^^^ 
is larger or smaller than the sound speed. In practice 
though, since both the sound speed and the viscosity 
damp the dark energy perturbations in a similar way, we 
can speak of only one effective sound horizon^ which, once 
crossed, suppresses the growth of dark energy perturba- 
tions; this makes c^ and c^.^g degenerate when considering 
dark energy density perturbations. 

We then have used the equations for Q and E to study 
analytically the effect of dark energy perturbations on 
the dark matter power spectrum, on the growth factor 
and on the ISW effect. The changes in the matter power 
spectrum and in the growth index at late times are of 
the order of a few percent on scales larger than the ef- 
fective sound horizon of the dark energy. We also have 
found that the ISW effect is modified mainly by the time 
variation of the parameter E' while modifications in E 
itself do not have appreciable consequences. The param- 
eter E is a combination of the dark energy density con- 
trast (through Q) and the anisotropic stress (through ry) 
which both prevent perturbations from growing. Because 
of this, the two effects are difficult to capture separately 
as they both act similarly. To disentangle the two effects 
is possible in principle by using joint observations which 
measure Q and E, such as e.g. the galaxy power spec- 
trum together with the weak lensing maps. In a follow-up 
paper [33| , we evaluate the impact of combining different 
experiments in order to measure at the same time the 
sound speed and the anisotropic stress of dark energy. 
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Appendix A: Decaying modes for the anisotropic 

stress 

In this appendix we look carefully at the full analytic 
solution of Sigma and in particular at the decaying modes 
which we have previously neglected. We start by looking 
at modes above the sound horizon; combining Eqs. ()12|) 
and (j2ip we obtain a second order differential equation 
for the dark energy anisotropic stress 
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being Sq the initial dark matter density contrast that acts 
as a source term for dark energy perturbations. The full 
solution, after some mathematical manipulation, is: 



Let us look more carefully: the two terms neglected in 
the solution Eq. (|A2p are a combination of Sin and Cos 
which arc oscillating functions if a > 0. Given that c^j^ 
is assumed to be non-negative in this work, a is always 
positive as long as — 1 < w < 0. In case of "phantom" 
dark energy instead there may be a growing solution as 
a may become negative. One way out of this problem 
is to assume a negative viscosity term only for phantom 
model, as pointed out in [21|. Moreover, the oscillating 
terms in Eq. (|A2|) are also decaying modes because they 
are multiplied by the term l/(aa)^. 



Let us now consider modes below the sound horizon, 
i.e. cl ^ 0; combining Eqs. (HH), dUl) and dH]), the full 
first order differential equation for the anisotropic stress 
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where ci and C2 are constants of integrations depending 
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where wc also include the term 9H a [c]. — w)/k (that 
we previously neglected as it is a decaying mode) and for 
simplicity wc set 0o = iH^^rno^Q/'^ ■ The full solution 
then reads 
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where ci is the integration constant related to the initial 
conditions of the anisotropic stress a. All terms inside the 
curly brackets except for the first one are decaying modes; 
moreover they all manifest a strong dependence on the 
scale k, becoming even less important at small scales. 



The last term in Eq. (jASp depends on the exponent 
24cL(c?-l-w;) 



8cL + 3c2(l + u;) 



(A6) 



The former is positive if c^^^ < or if c^ < 1 + w; even 
if the last case is allowed in our calculations wc should 
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remember that Eq. (jA5|) was evaluated for modes below 
the sound horizon, henee for c^ substantially larger than 
zero; if we allow c^ to be close to zero, then the quantity 
i6\\ becomes, for w ~ —0.8 



3(1 + w) =0.6 



(A7) 



which is a positive number; however the last term in 
Eq. (jASP is multiplied by an extra 1/a^ term which en- 
sures us that this is still a decaying mode. As c^ becomes 
larger than 1 + w, the only possibility for the last term 



in Eq. (jASp to be a growing mode is Cyj^ < 0. Again, 
in order to stabilise the entire solution if c^j^. becomes 
negative one has to allow w < —1, i.e. phantom dark 
energy. 



Appendix B: Decaying modes for the density 
contrast 



Here we notice that to distinguish between sound 
speed an the type of anisotropic stress used in this 
paper one cannot use observables which depend 
only on S but one also needs to add observables 
depending on velocity perturbations. 

• Modes above the sound horizon Also in this 
case we can write the solutions (|24p and (P5)) in 
terms of c^j^j:; however, the formers have been eval- 
uated under the assumption c^ ^ so that the 
effective sound speed now depends only on the vis- 



cosity term c.^jj, (and w) 



Ce//('^?=0) 



31 



(C3) 



Then it turns out that Eqs. dMl) and ^^ have the 



exact same form as Eqs. 



with ci = 0. 



Here we study the solutions for the dark energy density 
contrast using the results found in the previous section 
and verify our assumptions. We consider only modes 
above the sound horizon because it is the only case where 
a differential equation for the density contrast is involved. 
The solution is 



<5(a) 



30o(l + w) 
27F4fi?„o(l 



1 



2ci 



8cl^,k^w • 64cLfc^a(l + 3u;) 
+ a'"'c2 + 



c2 
w) 



V^i/ov/f^(l + 6«^) 



(Bl) 



where ci and C2 are the constants of integration and 0o = 
3i/o^mo'5o/2. The second and third term in the square 
bracket are decaying modes and hence can be neglected; 
the second term goes like a~^^^ and it is also a decaying 
mode, whereas the last term in Eq. (jBip is a growing 
mode if w > which is not the case of dark energy 



Appendix C: Solutions with the effective sound 
horizon 



In this appendix we write the solutions for the dark en- 
ergy density contrast and velocity perturbation in terms 
of the effective sound speed c^y ^ 

w)]. 



c2 + 8(c2-H/[3(l + 



• Modes below the sound horizon. Eqs. p7)) and 
(HHI now read 



S 

V 



(1 + w) (t>o 
3(1 + w) (c2 - u;) 



(CI) 
(C2) 
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